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SUMMARY 

The  aeroelastician  needs  reliable,  efficient  methods  for  the 
calculation  of  unsteady  blade  forces  in  turbomachines.  The  validity 
of  such  theoretical  prediction  models  can  only  be  established 
if  researchers  apply  their  flutter  and  forced  vibration  predictions 
to  a  number  of  well  documented  experimental  test  cases. 

The  purpose  of  the  present  work  is  to  select  a  certain  number  of 
experimental  standard  configurations  to  be  used  as  test  cases  for 
theoretical  prediction  models,  as  well  as  for  other  experiments. 

A  standarized  nomenclature  and  reporting  format,  to  be  used  for 
facilitation  of  comparisons  in  the  present  project,  has  been 
defined.  The  reporting  format  allows  for  both  detailed  comparison 
of  the  experimental  and  theoretical  results,  as  well  as  for  a 
direct  judgment  of  the  flutter  tendencies  of  a  specific  turboma¬ 
chine  blading. 

A  large  number  (36)  unsteady  experimental  data  have  been  compiled, 
classified  and  evaluated.  Out  of  these  36  data,  9  have  been 
retained  as  possible  test  cases. 

Also  the  theoretical  part  of  the  project  has  received  a  large  inte¬ 
rest.  This  is  manifested  by  the  fact  that  already  25  researchers 
(from  the  academic  and  industrial  worlds)  have  announced  their 
intentions  to  "predict"  the  aeroe lastic behaviour  of  one  or 
several  of  the  experimental  test  cases. 

The  potential  for  a  very  positive  result  from  this  joint  research 
work  has  become  increasingly  larger,  as  many  more  institutions 
than  originally  previewed  are  interested  in  the  establishment 
of  the  validity  of  existing  and  future  theoretical  prediction 
methods  for  flutter  and  forced  vibrations. 

In  total,  31  persons  are  participating  in  the  experimental  part 
and  25  in  the  theoretical.  The  participants  are  very  evenly 
distributed  throughout  the  Western  countries  with  an  interest  in 
this  field,  namely  18  from  the  USA  (7+11),  24  from  Europe  (16+8) 
and  14  from  Japan  (8+6). 
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proposal  on  Jan. 9, 1981. 


EPE-'.AUEAV.E  >-BO=-TO  =E  ; -1  'HE=’.’l2Ur  A-  •  lUU" 

1 .  Introduction 

In  axial-flow  turbomachines  considerable  dynamic  blade  loads 
may  occur  as  a  result  of  the  interaction  between  the  unstea¬ 
dy  flow  and  the  blades.  The  trend  towards  ever  greater  mass 
flows  or  smaller  diameters  in  the  turbomachines  leads  to 
higher  flow  velocities  and  to  more  slender  blades.  It  is  there¬ 
fore  likely  that  aeroelastic  phenomena,  which  concerns  the 
motion  of  a  deformable  structure  in  an  airstream,  will  in¬ 
crease  ever  more  in  future  turboreactors  (fan  stage)  and  in¬ 
dustrial  turbines  (last  stage) . 

The  large  complications,  and  high  costs,  of  unsteady  flow  mea¬ 
surements  in  actual  turbomachines  makes  it  necessary  for  the 
aeroelastician  to  rely  on  cascade  experiments  and  theoretical 
prediction  methods  for  minimizing  blade  failures  due  to  aero¬ 
elastic  phenomena.  It  is  therefore  of  great  importance  to  va¬ 
lidate  the  accuracy  of  flutter  and  forced  vibration  predic¬ 
tions  and  compare  theoretical  results  with  cascade  tests  and 
trends  in  actual  turbomachines. 

Several  well  documented  unsteady  experimental  cascade  data 
exists  throughout  the  world,  as  well  as  many  different  promi¬ 
sing  calculation  methods  for  solving  the  problem  of  unsteady 
flow  in  two-dimensional  and  quasi  two-dimensional  cascades. 
However,  due  to  different  basic  assumptions  in  these  methods, 
as  well  as  many  different  ways  of  representing  the  obtained 
results,  no  real  effort  has  been  made  to  compare  theoretical 
methods  with  each  other.  Further,  the  validity  of  these  theo¬ 
retical  prediction  analysis  can  only,  since  hardly  any  exact 
solution  are  known,  be  verified  by  comparison  with  experiments 
This  is  very  seldom  done,  partly  because  of  the  reasons  men¬ 
tioned  above,  partly  as  well  documented  experimental  data  nor¬ 
mally  are  of  proprietary  nature. 

It  is  the  purpose  of  the  present  project 
to  partly  remedy  the  above  mentioned  situation. 

This  can  be  done  by: 

(I)  establish  a  data  bank  of  test  cases  for  off-design  cal¬ 
culations 
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(II)  initiate  a  workshop,  in  which  specialists  can  compare 
and  discuss  their  results 

(III)  define  precise  reporting  formats  for  representing  the 
experimental  and  theoretical  results 

(IV)  initiate  additional  measurements  on  existing  test  faci¬ 
lities  . 

The  final  report  of  the  present  work  will  be  presented  at  the 

"Third  International  Symposium  on  Aeroelasticity  in  Turboma¬ 
chines",  to  be  held  in  Cambridge,  UK,  1984  (/ 9/ ,  /10/) 


< 
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As  described  in  detail  in  the  fist  semi-annual  progress  report 
(|7/),  most  of  the  major  institutions  in  the  field  of  aeroelasti- 
city,  both  in  USA,  Europe  and  Japan,  are  participating,  in  one 
way  or  another,  in  the  present  joint  research  work  (see  table  1). 
Out  of  the  34  sets  of  experimental  data  presented  in  table  1, 
the  ones  from 


ONERA  (Girault) 

ONERA  (Szechenyi) 

TOKYO  UNIVERSITY  (Tanaka) 

DETROID  DIESEL  (Jay) 

ALLISON 

UNITED  TECHNOLOGIES  (Carta) 

NASA  LEWIS  (Boldman/Ball) 

EPF- LAUSANNE  (Fransson/Sch ldf 1 i ) 

EPF-LAUSANNE  (Fransson/Schlhf li) 

NATIONAL  AEROSPACE  (Kobayashi) 

LAB.  TOKYO 


are  still  retained  as  possible  test  cases.  The  others  have  been 
eliminated  due  to  one,  or  several,  of  the  following  criteria: 

•In  order  to  exploit  to  the  maximum  the  comparison  between  theories 
and  experiments,  only  a  limited  number  of  test  cases,  concerning 
fully  aeroelastic  phenomena  in  turbomachine-cascades,  can  be 
choosen. 

Eventually,  rotor-stator  interactions  may  be  introduced  in  the 
present  work  at  a  later  time. 

•The  state-of-art  of  calculation  methods  must  be  taken  account  of. 
Therefore,  the  scope  of  the  present  joint  work  will  not  include 
fully  three  dimensional  experiments  and  calculations . 

•The  data  (at  least  some  preliminary)  must  be  ready  for  distribution 
towards  the  end  of  1982.  The  data  from  Whittle  Laboratory,  TH- 
Aachen  and  Toshiba-Kawasaki  will  unfortunately  not  fulfill  this 
criteria . 


>wf  l  dot*  w>d  gogg  »x*ct  joluUMfli  Iron,  which 


•  In  the  experiments,  several  parameters  should  be  varied. 

This  is  not  the  case  in  a  few  tests  directed  primarily 
towards  the  establishment  of  flutter  boundaries,  although 
these  tests  may  be,  as  in  the  case  of  the  ones  by  Brown 
Boveri  and  Westinghouse ,  very  well  documented. 

•  The  test  cases  should  treat  results  from  linear  or  annular 
cascades ,  not  real  machines  or  isolated  airfoils,  in  order 
not  to  diverge  the  scope  of  the  work. 

One  or  two  experimental  test  cases  will  be  designed  in  each  of 
the  four  main  velocity  ranges:  incompressible,  subsonic,  transonic 
and  supersonic.  The  compilation  of  these  data  are  currently  under 
way,  and  a  report  with  this  information  (|8{)  will  be  distributed 
to  the  participants  for  off-design  calculations  towards  the  end 
of  1982. 

2 . 2  Definition  of  precise  reporting  formats  for  the 
experimental  and  theoretical  results 

In  the  above  mentioned  compilation  of  experimental  data  as  test 
cases  for  off-design  calculations  (|8|),  the  first  chapters 
treat  the  standardization  of  the  aerodynamic  nomenclature  to  be 
used  by  all  participants  in  the  present  work.  This  is  a  very 
important  factor,  in  order  to  facilitate  the  evaluation  and 
validation  of,  as  well  as  the  comparison  between,  the  experimental 
and  theoretical  results. 

A  minimal  number  of  prescriptions  should  of  course  be  used  in  a 
joint  work  of  the  present  type.  However,  the  main  nomenclature 
must  be  standarized  in  order  to  avoid  misunderstandings  and  to 
facilitate  the  comparison  of  relevant  information.  The  notations 
recommended  are  based  upon  the  ones  in  references  | 1 1 1 o  |6|,  and 
especially  on  the  work  by  Carta  (  |1|,  |5|  ). 

2.2.1  Steady  Two-Dimensional  Cascade  Nomenclature 

The  profiles  under  investigation  are  arranged  ,  in  a  two-dimensional 
section  of  the  cascade,  as  in  Fig.  1.  In  this  figure,  all  the 
physical  lengths  are  scaled  with  the  chordlength  "c",  and  the 
nomenclature  in  table  2  is  used. 
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SYMBOL 

c 


s 

V 


P 


Subscripts 

6- 

i 

n. 

x 

? 

4 

Z 


♦  oo 


EXPLANATION 


DIMENSION 


chord  length  m 

time  averaged  pressure  coefficient  - 

incidence  angle  deg 

Mach  number  _ 

unity  vector  normal  to  blade  surface,  positive  inwards  _ 

2 

pressure  N/m 

dimensionless  vector  from  mean  pivot  axis  to  an  arbitrary 
point  on  the  mean  blade  surface  — 

unity  vector  tangent  to  blade  surface,  positive  in  positive 
coordinate  directions  — 

velocity  m/sec 

dimensionless  chord-wise  coordinate  — 

dimensionless  normal- to-chord  coordinate  — 


flow  angle  deg 
chordal  stagger  angle  deg 
dimensionless  blade  pitch  — 


center  of  gravity 

total  head  value 

pitch  axis  (see  Fig.  1) 

component  in  x-direction 

component  in  y-direction 

measuring  station  upstream  of  cascade 

measuring  station  downstream  of  cascade 

values  at  "infinity"  upstream 

values  at  "infinity"  downstream 


Superscripts 

(mM  mth  blade,  m=0,  t  1,1  2, _ If  the  amplitudes,  interblade  phase  angle,.. 

are  constant  for  the  blades  under  conside¬ 
ration,  this  superscript  will  not  be  used. 


/w 

(55) 


pressure  surface  of  profile 

steady  (=timc  averaged)  values.  This  superscript  will  only  be  used 

in  ambiguous  context. 

suction  surface  of  profile 
amplitude  of  complex  value 


TABLE  2  Steady  two-dimensional  cascade  nomenclature 
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2.2.2  Unsteady  Two-Dimensional  Cascade  Nomenclature 
Blade  Motion 


Fig.  2  is  a  schematic  representation  of  cascaded  two- 

dimensional  airfoils;  the  form  of  the  section  is  considered 

to  remain  rigidly  fixed  during  heaving  and/or  pitching 

oscillations(h* (x,y,t)  and  oc  (t)  ,  resp.  ).  The  components 

h  ,  h  and*  of  the  motion  vectors  h  and  ^  are  noted  in 
x  y 

complex  form  to  account  for  phase  differences  between  the  trans¬ 
lations  in  x  and  y-directions  as  well  as  the  rotation*. 


heaving  : 


— *  <“*)  - 4  ^ 

h  (x,y,t)=  h  (x,y)  e  (la) 


pitching 


It  is  here  assumed  that  the  torsional  motion,  for  the 
(m) th  blade,  preceeds  the  bending  motion  by  a  phase  angle 
Further,  if  the  vibration  frequency  is  identical  for 
all  blades,  the  superscript  (m)  on  will  be  dropped. 
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SYMBOL 

EXPLANATION 

DIMENSION 

A 

amplitude  At  k  for  pure  sinusoidal  heaving 

_  ^  _ 

-  **•  "  "  "  pitching 

deg 

-  or  Vdeg 

r;<o 

unsteady  force  coefficient  vector  per  unit  amplitude, 
positive  in  positive  coordinate  directions 

clCt) 

unsteady  lift  coefficient  per  unit  amplitude 
positive  in  positive  y-direction 

-  or  Vdeg 

unsteady  momemnt  coefficient  per  unit  amplitude, - 
positive  in  clockwise  d'  cion 

-  or  /deg 

or  Vdeg 

V‘,0 

unsteady  pressure  coefficient  per  unit  ■amplitude 

C  gy 

aerodynamic  work  done  on  the  system  during  one  cycle 
of  oscillation 

-  or  Vdeg 

f 

vibration  frequency 

Hz 

td) 

dimensionless  (with  chord)  bending  vibration, 
positive  in  positive  coordinate  directions 

- 

~k 

dimensionless  bending  amplitude 

Tc 

reduced  frequency  = 

S  i-r 

t 

oUA) 

At, 


*  ' 


e 


ft',*) 


<t>r 

<Pu 


4>fi*) 

tA) 


(? 


Strouhal  number 
time 

pitching  vibration,  positive  nose  up 

pitching  amplitude 

pressure  difference  coefficient 

phase  lead  of  pitching  motion  towards  heaving  motion  of 
blade  (m) 

aerodynamic  damping  coefficient,  positive  for  stable  motion 

phase  lead  of  force  coefficient  towards  motion 

phase  lead  of  lift  coefficient  towards  motion 

phase  lead  of  moment  coefficient  towards  motion 

phase  lead  of  pressure  coefficient  towards  motion 

phase  lead  of  pressure  difference  coefficient  towards  motion 

circular  frequency  =  <2  Tff- 

intcrblade  phase  angle  between  blade  "CT  and  blande  "m". 

for  constant  interblade  phase  angle.  d*"*Jis  positive 


when  blade  "m"  leads  blade  "0" 


sec 

deg 

deg 


1, 


or  /deg 
deg  or  rad 


-  or  /deg 
deg  or  rad 
deg  or  rad 
deg  or  rad 
deg  or  rad 


rad/ sec 
deg  or  rad 


Subscripts 

A 

A*  h  for  heaving 

Superscript 

(pj) 

for  pressure  surface 

A*i  for  pitching 

(*») 

for  suction  surface 

I 

imaginary  part 

(m) 

blade  number 

R 

real  part 

m=  0,  £  l,t  2,... 

TABU:  3 


Unsteady  two-dimensional  cascade  nomenclature 
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Two-Dimensional  Aerodynamic  Coefficients 

The  unsteady  (complex)  blade  surface  pressure  coefficient, 

as  well  as  the  lift,  force  and  moment  coefficients,  are  scaled 

with  the  amplitude  of  the  corresponding  motion  (amplitude = A, 

-  -<-l 

where  A  =  h  or  <*.  ) . 


<«> 

«  cr  (< i) 

rA 


i .  f  I 

*  -  y>~  J 


U) 


y>U,i)  = 

profile 

Su.T*f*k.^C 


where 


(S’) 


profile 

Sur#«» 


force  components  are  positive  when  acting  in  positive  coor¬ 

dinate  directions. 

‘  is  positive  when  acting  in  clockwise  direction 
(B)  superscript)  denotes  the  blade  pressure  surface  (ps) 

^  or  blade  suction  surface  (ss). 


A  further  important  quantity,  for  slender  blades,  is  the  dimen¬ 
sionless  pressure  difference  on  the  blade,  ACp  .  This  is  defined 
as  the  difference  between  the  pressures  on  the  blade  pressure 
and  suction  surfaces; 
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AC|/X»4)  '  {Cfi  <X/0  ~ 


(6) 


All  of  the  above  mentioned  variables  can  be  expressed  in 
either  complex  exponential  form  or  in  component  form  as: 
.  *  S  (7a) 

w  *  I 

t  tVX 


-{■ 


Pw. 


£*)  *  '<■  t* 


)} 


(7b) 


Here,  the  subscripts  "R"  and  "I”  denotes  the  real  and 
imaginary  parts  of  the  pressure  coefficient  c  .  Physically, 
these  two  parts  can  be  interpreted  as  the  components  of  the 
pressure  coefficient  which  are  in-phase  (real  part)  and  out-of¬ 
phase  (imaginary  part)  with  the  blade  vibration.  Further,  the 
phase  angles  4,  are  all  defined  positive 

when  the  pressure  (pressure  difference,  lift,  force  or  moment, 
resp.)  leads  the  motion. 

The  amplitude  and  phase  relationships  in  eq.  (7)  are  defined  in 
the  usual  manner,  that  is: 
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It  should  here  be  noted  that,  in  computing  the  blade  surface 
pressure  distribution,  only  components,  and  not  amplitudes  or 
phase  angles  may  be  differentiated  (|1|).  Therefore 
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Two-Dimensional  Aerodynamic  Work 


The  two  -dimensional  differential  work  done  on  a  rigid  system  b^ 
the  aerodynamic  forces  and  moments  is  conventionally  expressed  by 
the  product  of  the  real  parts  (  in  phase  with  motion  components) 
of  force  and  differential  translation,  as  well  as  moment  and 
differential  torsion.  Thus,  the  total  aerodynamic  work  coefficient 
per  period  of  oscillation,  done  on  the  system  is  obtained  by  com¬ 
puting 


c,,  = 
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where 
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In  the  case  of  pure  sinusoidal  normal-to-chord  bending  or 
torsional  vibration,  as  well  as  sinusoidal  lift  and  moment 
responses,  respectively,  the  expressions  (11)  may  be  integrated 
into  the  following  simple  formulas 


cu/K  ~  ^  ^  c/_i  (r  c l  <^u  1 

Coc/flc  r  (=  O 


The  airfoil  damps  therefore  the  motion  when  the  imaginary  part 
of  the  lift  and  moment  coefficients,  resp.  is  negative. 

An  aerodynamic  damping  parameter  defined  in  references 

111,  as, with  the  same  assumption  as  in  (11), 


This  parameter  is  therefore  positive  for  a  stable  motion. 


All  theoretical  prediction  methods  for  flutter  and  forced 
vibrations  available  today  make  a  few  basic  assumptions. 

Most  of  these  methods  are  submitted  to  restrictions  regarding 

•  sinusoidal  blade  vibrations 

•  sinusoidal  pressure  responses 

•  identical  vibration  frequencies  for  all  blades 

•  identical  vibration  amplitudes  for  all  blades 

•  constant  interblade  phase  angles 

In  experiments,  however,  these  assumptions  can  never  be  cxactely  fullfilled. 
The  large  energy  hput  needed  to  drive  a  cascade  with  prescribed 
frequencies,  amplitudes  and  phase  angles  makes  it  impossible  to 
satisfy  the  three  latter  assumptions,  apart  from  tests  with 
low  frequencies  and/or  small  amplitudes.  Even  in  this  case  though 
the  pressure  response  on  the  profiles  will,  in  generally,  not  be 
sinusoidal . 

In  the  computation  of  the  aerodynamic  work  and  damping  coefficients, 
this  non-sinusoidal  pressure  response  on  the  vibrating  blades 
can  still  be  used  relatively  easy,  as  it  can  be  shown  that  only 
the  frequency  of  the  pressure  response  spectra  that  corresponds 
with  the  blade  vibration  frequency  contribute  to  the  aerodynamic 
work . 

However,  non-identical  blade  vibration  amplitudes  and  interblade 
phase  angles  may  largely  contribute  to  some  discrepancies  between 
the  experimental  and  theoretical  results, 

A  detailed  comparison  and  evaluation  of  unsteady  experimental  and 
theoretical  cascade  results  necessitate  thus  a  knowledge  of  the  ac¬ 
tual  time  histories  of  the  experiments. 

From  these  time-recordings,  a  statistical  evaluation,  or  a 
Fourier  analysis  can  be  used  to  appreciate  how  well  the  different 

idealizations  in  the  prediction  models  approximate  the  real 
physical  conditions. 

In  the  final  comparison  between  the  experimental  and  theoretical 
results,  such  an  evaluation  will  therefore,  if  it  exists  for  the 
experiment  under  consideration,  be  presented. 
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2.2.3  Definition  of  precise  reporting  formats 

It  is  obvious  that  no  universal  representation  of  the  flutter 
phenomena  can  be  defined  as  long  as  the  underlying  physical 
reasons  for  the  self-excited  vibrations  are  not  understood. 
Secondly,  different  importance  is  attached  to  different  purposes 
of  the  experimental  and  theoretical  investigations.  It  is 
evident  that  a  designer  of  a  turboreactor  or  industrial  turboma¬ 
chine  is  not  interested  in  the  same  details  of  the  aeroelastic 
phenomena  as  an  experimentalist,  who  is  performing  tests  in, 
for  example,  an  annular  cascade.  Further,  a  numerical  analyst 
has  another  sort  of  problem  (mesh  generation,  accuracy...)  and 
tend  therefore  to  present  his  results  in  a  way  to  facilitate 
comparisons  with  other  theoretical  approaches. 

One  of  the  main  objectives  of  the  investigations  for  both 
experimentalists  and  theoreticians  is,  of  course,  to  furnish 
the  designer  of  turbomachines  with  enough  information  to  enable 
him  to  minimize  blade  failures  due  to  aeroelastic  phenomena. 
However,  such  a  universal  method  is  still  years  to  come. 

The  representations  used  in  the  present  comparative  work  will 
therefore  allow  for  evaluations  and  illustrations  between 

•  experimental  -  experimental  results 

•  experimental  -  theoretical  results 

•  theoretical  -  theoretical  results 

It  will  also  give  a  possibility  of  direct  comparison  of  the  flutter 
tendencies  of  the  different  cascades  under  investigation. 

Thus,  two  different  reporting  formats  of  the  results  will  be  used: 

I  The  first  part  is  directed  towards  the  comparison 
and  establishment  of  validity  of  theoretical  models 
(and  the  experimental  result).  It  will  also  indicate 

a  quantitative  comparison  between  the  steady  and  unsteady 
blade  forces. 

II  The  second  part  concerns  the  direct  use  of  the  data  for 
establishing  of  the  flutter  boundaries. 


* 
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In  both  cases,  the  exact  representation  will  largely  depend 
upon  the  existing  experimental  data. 

Further  to  these  representations,  every  participant  is  asked  to 
present  his  data  in  any  other  way  he  judges  useful  for  better 
physical  understanding  of  the  flutter  phenomena. 

I  Detailed  comparison  of  experimental  results  and  theoretical 
approaches 

The  establishment  of  the  validity  of  both  experimental  and 
theoretical  results  can  only  be  done  by  a  mutual  agreement 
between  the  measured  and  calculated  blade  pressure  distributions. 
This  detailed  comparison  between  experiments  and  theories  will 
therefore  be  done  on  the  basis  of  diagram  1,  which  will  be  presented 
for  a  few  different  combinations  of 

•  interblade  phase  angle 

•  reduced  frequency 

•  inlet  conditions 

•  cascade  geometry 

depending  upon  the  existing  experimental  data  for  the  test  case 
under  investigation. 

Quite  a  few  prediction  models  for  flutter  or  forepd  vibrations 
are  based  upon  small  perturbation  theories,  which  means  that  the 
steady  pressure  distribution  on  the  blade  must  be  given  as  an  input 
data.  In  this  diagram  1,  the  time-averaged  blade  pressure  distribution 
will  therefore  be  specified  for  such  studies. 

Further  the  comparison  between  the  steady  and  unsteady  blade  pressu 
re  distributions  may  in  some  cases  give  a  quantitative  notion 
about  the  aeroelastic  phenomena  under  investigation  (instabilities 
due  to  stall,  choke,  shock...  ,  coupling  effects  between  the 
steady  and  unsteady  flow  fields). 

The  difference  of  the  pressure  coefficients,  (*)  ,  along 

the  blade  indicates  the  distribution  of  the  stable  and  unstable 
parts  of  the  blades.  This  coefficient  is  therefore  also  of  interest, 
and  will  be  plotted  in  diagram  2. 


DIAGRAM  |  A  t  STEADY  BLADE  PRESSURE  COEFFICIENT 


DIAGRAM  2  i 


MAGNITUDE  ANO  PHASE  LEAD  OF  PRESSURE 
DIFFERENCE  COEFFICIENT  ON  THE  BLADE 
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1 1 .  Flutter  boundaries 

This  second  illustration  concerns  the  integrated  values  of 
the  aerodynamic  blade  forces  and  moments,  as  well  as  the 
aerodynamic  work  and  damping  coefficients. 

Two  different  representations  (see  diagram  3  and  4)  will  be 
used  to  elaborate  the  flutter  dependancies  of  several  important 
parameters 

•  reduced  frequency 

•  interblade  phase  angle 

•  inlet  flow  velocity 

•  inlet  flow  angle 

•  outlet  flow  velocity 

•  cascade  geometry 

First,  the  unsteady  blade  pressure  coefficients  will  be  integrated 
to  yield  the  aerodynamic  lift  and  moment  coefficients  as  in 
diagram  3.  The  phase  angles  <j>^  and  res p. ,  give  then  directly 
information  about  the  aeroelastic  stability  of  the  system 
(see  chapter  2.2.2) . 

Secondly,  the  aerodynamic  work  anddimping  coeeficienti  per 
cycle  of  oscillation  may  be  calculated  if  the  mode-shape  of  the 
motion  is  well  known.  Most  of  the  problems  treated  in  the  present 
work  will  concern  rigid  systems  (at  least  for  the  theoretical 
predictions),  wherefore  the  aerodynamic  damping  coefficient  can 
be  easily  computed  and  plotted.  This  coefficient  is  the  single 
most  interesting  in  judging  the  aeroelastic  behaviour  of  a 
specific  cascade,  and  should  therefore  be  plotted  for  the  benefit 
of  the  turbomachine  designer. 
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2 . 3  Theoretical  off-design  calculations 


It  is  obvious  that,  for  a  detailed  evaluation  of  the  state-of-art 
of  prediction  models  for  aeroelasticity ,  several  codes  must 
be  compared.  To  this  end,  52  institutions  have  been  contacted 
regarding  the  off-design  calculations  of  the  test  cases. 

Up  to  now  the  following  25  persons  have  announced  their  partici¬ 
pation,  each  with  one  or  several  computing  codes  (see  table  4). 

The  off-design  calculations  will  be  performed,  as  suggested  by 
the  Scientific  Committee  at  the  Lausanne  Symposium  on  Aeroelastici¬ 
ty  in  Turbomachines  J 1 0 1  ,  by  the  participants  during  1983  ,  where¬ 

after  all  the  results  will  be  compiled  and  compared  for  presenta¬ 
tion  at  the  third  Symposium  on  Aeroelasticity,  to  be  held  in 
Cambridge,  U.K.,  in  1984. 


Institution 


Name 


Number  of 
Codes 


(Physical  Sciences  Inc.} 
^University  of  Tokyo  J 

USA 

Japan 

N.H.  Kcmpl 

H .  Sho j i  j 

1 

Cambridge  University 

UK 

D.S.  Whitehead 

4 

United  Technologies 
Research  Center 

USA 

J.M.  Verdon/ 

J.M.  Caspar 

1 

United  Technologies 
Research  Center 

USA 

J.M.  Verdon 

1 

University  of  Tokyo 

J  apan 

S .  Ka j i 

1 

ON  ERA 

France 

P.  Salaiin 

3 

NASA  Lewis  Research 
Center 

USA 

M.F..  Goldstein/ 
W.ll.  Braun/ 

F.B.  Molls 

2 

Kuyushu  University 

J  apan 

M.  Namba 

1  (2?) 

Florence  University 

I  taly 

F.  Martel li 

2 

Lausanne  Inst,  of 
Technology 

Switzerland 

T.H.  Fransson 

1 

University  of 

Notre-Dame 

USA 

H.  Atassi 

9 

DFVLR-G5tt ingen 

W.  Germany 

V.  Carstens 

9 

University  of  Tennes¬ 
see  Space  Inst. 

USA 

J.  Caruthers 

9 

Technische  Hochschule 
Aachen 

W.  Germany 

H.E.  Gallus/ 

K.  Vogeler 

9 

Massachusetts  Inst, 
of  Technology 

USA 

E.  Crawley 

3 

Imperial  College 

UK 

M.  Graham 

? 

University  of  Tennessee 
Space  Institute 

USA 

M.  Kurosaka 

9 

University  of  Tokyo 

Japan 

H.  Tanaka 

9 

University  of  Tokyo 

Japan 

Y.  Tanida 

9 

Naval  Postgraduate 
School 

USA 

M.  Platzer 

9 

Mitsubishi  Heavy 

Ind.  Ltd. 

Japan 

S.  Takahara 

9 

Table  4:  Participants  in  the  theoretical  part  of  the  work 

on  "Aeroe 1  as t ic  i  ty  in  Turbomach i ne-Cascadcs"  * 
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The  observations  during  the  present  project  have  shown  that 
experimental  aeroelastic  investigations  in  highly  loaded 
subsonic  turbine  bladings  in  plunging  and  pitching  modes  are  diffi¬ 
cult  to  find  in  the  open  literature.  To  remedy  this,  the  Labora- 
)ire  de  Thermique  AppliquSe  (LTA)  has  started  such  an  experimental 
work.  The  cascade  is  currently  under  fabrication,  and  it  will  be 
instrumented  with  flush  mounted  pressure  transducers  and  static 
pressure  tappings  towards  the  end  of  this  year.  It  is  expected 
that  the  first  aeroelastic  results  with  this  cascade  will  be 
obtained  in  March-April  1983,  whereafter  these  data  will  be  dis¬ 
tributed  as  test  cases  for  the  off  design  calculations. 

Another  lack  of  well  documented  experimental  data  have  been  found 
in  the  supersonic  flow  regime.  It  is  here  well  known  that  the 
linearized  supersonic  flat  plate  theories  gives  very  good  agree¬ 
ments  with  two-dimensional  compressor  cascade  experiments  (thin 
airfoils).  However,  the  validity  of  the  strip  theory,  i.e.  the  use 
of  two-dimensional  prediction  methods  in  a  quasi-two-dimensional 
or  a  three-dimensional  upstream  flow,  has  never  been  established 
in  the  transonic  and  supersonic  flow  regimes,  not  even  for  flat 
plates.  Theoretical  investigations  by  Namba  |ll|  indicate  that 
three-dimensional  effects  in  completely  supersonic  cascades  are 
generally  small  and  the  strip  theory  predicts  local  as  well  as 
integrated  aerodynamic  forces  with  good  accuracy.  In  the  transonic 
flow,  i.e.  in  the  normal  operating  range  of  a  modern  turboreactor, 
however,  the  strip  theory  approximation  breaks  down  near  the 
sonic  span  station  and  three-dimensional  effects  are  of  primary 
importance.  An  experimental  investigation  of  this  phenomena  in 
a  quasi-two-dimensional  test  facility  would  therefore  be  of  large 
interest  for  the  development  of  fan  stages  in  future  jet  engines. 
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3.  RESEARCH  PLANNED  FOR  THE  PERIOD  SEPTEMBER  1,  1982  - 
AUGUST  31,  1983 

As  can  be  seen  from  chapter  2,  the  present  project  is  now  well 
under  way. 

During  the  rest  of  1982,  refinements  in  the  compilation  of  test 
cases  will  be  made.  These  will  thereafter  be  distributed  for  off- 
design  calculations  by  several  prediction  methods  towards  the 
end  of  1982. 

The  theoreticians  participating  in  the  present  work  will  thereafter 
"predict"  the  aeroelasticity  of  the  test  cascades.  These  calcula¬ 
tions  will  all  be  performed  during  1983. 

Further,  many  contacts  will  be  necessary  in  order  to  supply  the 
theoreticians  with  all  the  information  each  one  need  for  their 
prediction  methods,  as  well  as  to  evaluate  and  compare  the 
different  results  obtained  from  the  different  experiments  and  theo¬ 
retical  models. 
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4.  FISCAL  STATUS  FOR  THE  FIRST  YEAR  OF  THE  PROJECT 
September  1,  1981  -  August  31,  1982 


This  chapter  is  divided  into  two  parts: 

I  Coordination  of  the  present  joint  research  project 

II  Contributions  to  aeroelast ici ty  research  by  LTA 

(See  Unsolicitated  Research  Proposal  -  Jan.  1981,  chapter 
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4 . 1  Coordination  of  the  present  joint  research  project 


AFOSR  Costsharing 


actual 
(12  months) 

US$ 


budgeted 
(12  months) 

US$ 


actual 
(12  months) 

US$ 


budgeted 
(12  months) 

US$ 


1.  Salaries  &  Wages 

-Principle  Investigator 
P.  Suter 


-Research  associate 
T.  Fransson  11' 200. — 


2.  Expendible  Material 
&  Supplies 

3.  Permanent  Equipment 

4.  Travels 


1  trip  to  Germany 
(KWU,  DFVLR.Th- Aachen) 


1  trip  to  G.B. 
(ASME-Conference , 
Rolls-Royce,  GEC, 
Central  Elec.  Gen. Board, 
Cambridge  University) 


1 '600. — 


1  trip  to  France 
(ONERA) 

5.  Publication,  report, 
initialization  of  the 
project 

Secretary,  telex, 
telephone 

6.  Other  direct  costs 


7.  Indirect  costs 


7 '000. — 

9 ' 700. —  6 ' 500 . — 


3 '  100 .  — 

480. — 
300.— 

220.— 


2 ' 200 . — 
3 ' 500 . — 


6 . 600 . — 


5'000. — 


1 ' 800.— 


l'OOO.— 

S'000.— 


Total  costs  12'800. —  12'800.—  20’200.— 


19'  400.— 
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4.2  Experimental  and  theoretical  acroelastic  investigations 
at  LTA 

Approximative  1982  budget  of  turbomachine  section  at  LTA  (for 
information  only) 


1. 

Salaries  &  Wages 

7  collaborators,  about  3,9  man/year 

1 10 '000. 

2. 

Expendable  materials 

68'000. 

3. 

Computer  time 

42*000. 

4. 

Other  direct  costs 

8*000. 

Approximative  1982  budget  of  the  turbomachine 

section  at  LTA  $  228'000. 


5.  CONCLUSIONS 


In  the  first  year  of  the  project  "Aeroelasticity  in  Turbomachine- 
Cascades",  a  very  large  interest  for  the  present  work  have  been 
found  to  exist  both  in  the  academic  world  as  among  the  manufacturers 
of  turboreactors  and  industrial  turbomachines.  Unfortunately,  several 
industrial  institutions  cannot  participate  because  of  the  proprie¬ 
tary  nature  of  their  work.  Still,  a  large  number  of  experimental 
data  (36  sets)  have  been  put  at  our  disposal  for  comparison  between 
experimental  results  and  prediction  models  for  flutter  and  forced 
vibrations . 

Also  the  theoretical  aspects  of  the  comparative  work  is  well  represen¬ 
ted,  as  already  2S  persons  have  announced  their  intentions  to 
participate  with  one  or  several  prediction  methods. 

From  the  large  date  base  of  experimental  results,  a  few  test  cases 
are  currently  being  compiled  for  distribution  for  off-design 
calculations,  and  comparisons,  towards  the  end  of  1982. 

As  a  part  of  this  compilation,  a  precise  nomenclature  and  representa¬ 
tion  format  to  be  used  during  the  present  work  has  been  established. 

An  absence  of  well  documented  experimental  data  have  been  found, 
noticably  in  the  domain  of  high  turning  subsonic  turbine  blades 
and  in  the  field  of  transonic/supersonic  quasi-two-dimensional 
flat  plate  investigations.  As  regards  to  the  high  turning  turbine 
blades,  the  Laboratoire  de  Thermique  Appliquee  has  reocntlv  started 
such  experimental  work.  The  result  of  this  investigation  can 

be  distributed,  as  a  further  test  case,  during  1983. 
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A.  Ishikawa  of  Oscillating  Supersonic  and  Transonic  Annular 
Cascades" 

ASME-Paper  82-GT-126  ,  1982 

D.  Schldfli  "Flatterunte rsuchungen  an  einem  t ranssonischen 

Turbinengitter  im  Ringgitterkanal  an  der 
EPF-Lausanne" 

Report  EPFL-LTA/TM-5-82 
Lausanne ,  Switzerland ,  1982 

W.J.  Rakowski  "A  Research  Program  for  the  Experimental 
Analysis  of  Blade  Stability" 

AIAA/SAE  -  14  joint  Propulsion  Conference , 

Las  Vegas ,  1978 

M.  Kurosaka  "Linear  and  Nonlinear  Analysis  of  Vortex 
Whistle-Another  Blade  Buster” 

in  1 10\  ,  pages  443-453 

A.  Kirschner  "Control  of  Vibration  in  Aeroelastic  Cascade 

B.  Fosco  Experiments" 

E.  Muller  |i0)  ,  pages  285-295 

A.  Boles  "A  Test-Facility  for  the  Experimental  Inves¬ 

tigation  of  Steady  and  Unsteady  Flows  in  Annu¬ 
lar  Cascades" 

To  be  presented  at  the  ASME  Gas  Turbine  Con¬ 
ference ,  1983. 

T.H.  Fransson  "Numerical  Experiments  with  Different  Boundary 
Conditions  in  Subsonic  Unsteady  Flows" 

Presented  at  "Unsteady  Perturbations  in  Inter¬ 
nal  Flou8"  -  Cambridge ,  England, 

19-20/3  1981 

T.H.  Fransson  "Numerische  Berechnung  der  zweidimens ionalen 
Unterschallstromung  in  einem  os zi 1 1 ierenden 
Schaufelgitter" 

Presented  at  the  "Stromungstag  an  der  EPF- 
Lausanne",  Switzerland,  15.01.1982 


Re;  Aeroelastic  research  work  at  LTA  (Lab-  de  thermique  appliqufie) 

LTA  is  working  actively  in  the  field  of  aeroelasticity  with 

two  different  projects,  both  of  which  will  be  put  at  the  dis¬ 
posal  of  the  workshop: 

•  experimental  flutter  research. 

•  theoretical  forced  vibration  investigations. 

I . i .  Experimental  work 

•  The  experiments  with  an  annular  transonic  turbine  cascade 
(12,15,16)  have  permitted  us  to  observe  self-excited  blade 
vibrations  in  our  non-rotating  annular  test  facility.  This 
instability  can  be  classified  as  shock-  and/or  stall-induced 
flutter  (see  fig.  1) . 


Outlet  Mach  number 


Fig  1:  Flutter  boundaries  for  transonic  turbine  cascade  at  LTA 


It  has  been  found  that  flutter  is  an  organized  movement,  in 
which  the  amplitude  of  a  certain  blade  fluctuates  in  time. 

It  is  seen  that  a  continuous  energy  transport  from  the  first 
to  the  second  etc.  to  the  last  blade  in  the  cascade  takes  pla¬ 
ce.  This  coupling  over  all  blades  in  an  annular  nozzle  makes 
it  rather  unlikely  that  real  flutter  conditions  may  be  repro¬ 
duced  in  a  linear  cascade. 

At  the  moment,  this  transonic  turbine  cascade  is  run  under  forced 
vibrations.  Systematics  tests  are  performed  to  establish  the 
influence  of  several  important  parameters,  as  for  example  the 
interblade  phase  angle  and  backpressures,  on  the  aerodynamic 
stability  of  the  cascade. 

As  the  next  step  in  aeroelastic  experiments  at  LTA,a  highly 
loaded  subsonic  turbine  cascade  for  forced  vibration  experi¬ 
ments  is  under  fabrication  and  will  be  built  into  the  non¬ 
rotating  annular  test  facility.  The  objectives  are  threefold: 

+  gain  experience  with  a  highly  loaded  subsonic  cascade 
in  unsteady  motion 

+  result  for  off-design  calculations  (both  for  the  calculations 
performed  at  LTA  and  for  the  working  group  on  acroelast  icity 
in  general) 

+  compare  results  obtained  in  an  annular  cascade  with  correspon¬ 
ding  ones  in  a  linear  cascade.  To  this  end,  we  are  looking 
for  a  laboratory  equiped  with  a  linear  cascade  for  forced 
and  self-excited  vibrations,  in  order  to  coordinate  our  research 
and  perform  experiments  on  an  identical  cascade. 

Among  the  results  to  be  obtained  are: 

+  steady  and  unsteady  pressure  probe  traverses  up-  and  down 

+  steady  and  unsteady  pressure  coefficients  on  blades 

+  b lade  vibrat ion 

+  vizualization  in  the  transonic  flow  domain  (Schliercn) 

+  instability  dependance  on  flow  conditions,  b 1 ade- vibrat ions , 
interblade  phase  angle,... 


The  instrumentation  (unsteady  flushmountcd  pressure  transducers) 
will  follow  later  this  year,  and  the  first  tests  are  scheduled 
for  March  1983. 


Several  other  interesting  phenomena  have  been  investigated  in 
this  annular  test  vehicle.  As  an  example  can  be  mentioned 
that,  in  certain  flow  domains,  the  same  "Vortex  Whistles" 
as  was  noticed  in  the  annular  test  rig  at  General  Electric  Co, 
appeared  also  in  our  nozzle  (  Jl 3j  )  |l4|  ) 

In  our  case,  this  problem  could  be  remedied  simply  by  repro¬ 
ducing  the  same  flow  conditions  with  another  setting  of  the 
upstream  swirl  vanes. 

Experimental  investigations  on  blade  vibrations  due  to  ran¬ 
dom  fluctuations  of  the  upstream  flow  conditions  are  also 
under  way  at  LTA.  These  tests  are  performed  on  highly  loaded 
subsonic  turbine  blades  in  a  linear  cascade. 

Among  the  most  interesting  results  can  be  mentioned: 

+  Miniaturized  pressure  transducers  mounted  on  the  blades 
indicate  a  high  correlation  between  the  blade  movement 
and  the  pressure  variation  in  the  attached  flow  region. 

+  In  the  detached  flow  domain,  however,  the  pressure  fluctuations 
do  not  show  a  correlation  with  the  blade  movement. 

+  A  high  upstream  flow  turbulence  excites  the  blades  into 
a  much  larger  amplitude  than  a  lower  turbulence  level. 


* 
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1 . 1 1 .  Theoretical  work 

A  numerical  code  for  predicting  instabilities  in  a  subsonic 
cascade  with  vibrating  blades  is  being  developed  at  the  LTA. 

This  work  has  shown  that  the  inlet  and  outlet  boundary  conditions 
may  influence  the  unsteady  flow  in  the  cascade  to  a  large 
extent.  Some  results  of  this  calculations  have  been  represented 
in  a  movie  film,  showing  clearly  how  the  pressure  waves  created 
at  the  vibrating  blades  travel  up  -  and  downstream,  and  how  they 
interfere  with  the  upstream  and  downstream  boundaries  (  / 6 /  ,  /!/). 
The  code  is  currently  being  changed  to  allow  better  for  higher 
stagger  angles  and  thicker  blades.  This  improvement  is  judged 
necessary  in  order  to  correctly  handle  the  unsteady  flow  in 
the  leading  and  trailing  edge  regions. 
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I. Ill  Scientific  reports  since  the  submission  of  the  research 
proposal  on  January  9,  1981 


•  A.  Boles  "Berechnung  der  Stromung  in  einem  Ueberschall- 

angestromten  Schaufelgitter  bei  variabler 
Schichtdicke" 

Z AMP  Vo l .  32. 1981 


•  A.  Boles 


•  T.  Fransson 


•  T.  Fransson 


•  A.  Boles 
A.  Mathys 


•  D.  SchlSfli 


•  T.  Fransson 


•  P.  Testori 
A.  Schiess 


"Charakteristikenverfahren  fur  quasi- 
dreidimensionale  und  anisentrope  Ueber- 
schal lstromungen" 

Presented  at  "Stromung stag  an  der  EPFL" 
Lausanne ,  Switzerland,  IS. 01. 1982 

"Numerical  Experiments  with  Different  Boun¬ 
dary  Conditions  in  Subsonic  Unsteady  Flows" 

Presented  at  "Unsteady  Perturbations  in  Inter¬ 
nal  Flows" 

Cambridge,  England,  19/20.3.1981 

"Numerische  Berechnung  der  zweidimensionalcn 
Unterschal lstrbmung  in  einem  osz i 1 1 icrenden 
Schaufelgitter" 

Presented  at  " Stromung s tag  an  der  EPFL" 
Lausanne,  Switzerland,  IS. 01. 1982 

"Experience  with  Turbine  and  Compressor 
Cascades  in  an  Annular  Test  Facility" 

Proceedings  of  Symposium  on  "Measuring  Tech¬ 
niques  in  Transonic  and  Supersonic  Cascades 
and  Turbomachines” , 

Lyon,  France,  1981 

"Flatteruntersuchungen  an  einem  transsonischen 
Turbinengitter  im  Ringgitterkanal  an  der 
EPF-Lausanne" 

Report  EPFL-LTA/TM-S-82 ,  Lausanne,  Switzerland 
1982 

"Abloseverhalten  eines  Kreisbogcnprofi Jes 
bei  variabler  Reynoldszahl  und  variablem 
Anstrdmwinkel" 

Report  EPFL-LTA/TM- 1-82 

"Etude  d'^coulement  autour  d'une  aubc  isol6e 
en  r6gime  transsonique  (Diplome  1981)" 


Report  EPFL-LTA,  1981 


F.  Beretta-Piccoli 
A.  BSlcs 


"Versuche  mit  einem  Radial vcrdichter- 
Vorschaltrad:  2.  Schlussbericht" 

Report  EPFL-LT A/TM-6-82 


4  T.  Fransson 
0.  Sari 


•  F.  Beretta-Piccoli 
T.  Fransson 


•  X.C.  Lian 


•  0.  Sari 

F.  Beretta-Piccoli 
T.  Fransson 

•  T.  Fransson 
P.  Suter 


•  T.  Fransson 
X.  C.  Lian 


•  T.  Fransson 
P.  Suter 


"Charasteristics  of  Aerodynamic  Five- 
Hole-Probes  in  Transonic  and  Supersonic 
Flow  Regimes” 

Proceedings  of  "Measuring  Techniques 
in  Transonic  and  Supersonic  Cascades  and 
Turbomachines " 

Lyon ,  France ,  15/16.10.81 

"Probe  Investigations  in  the  Proximity 
of  a  Wall  in  Supersonic  Flow" 

Proceedings  of  "Measuring  Techniques 
in  Transonic  and  Supersonic  Cascades  and 
Turbomachines " 

Lyon,  France,  15/16.10.81 

"Visualisation  de  l’ecoulcment  de  surface 
sur  un  profil” 

Report  EPFL-LTA/TM-2-82 

"Hichung  der  Aerodynamischen  Sonde  R  18" 
Report  EPFL-LTA/TM-Z-82 


"Aeroelasticity  in  Turbomachine-Cascades" 
First  Semi-Annual  Progress  Report 

Report  EPFL-LT A/TM-4-82 


"Ablflseverhalten  eines  Kreisbogenprofiles 
bei  variabler  Reynoldszahl  und  variablcm 
Anstrftmwinkel" 

(Nachtrag  zum  Bericht  LTA-TM-1-82) 

Report  EPFL-LTA/TM-7-82 


"Aeroelasticity  in  Turbomachine-Cascades" 
Second  Semi-Annual  Progress  Report 

Report  EPFL-LT A /TM- 8 -82 
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